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PROCEDE ET DISPOSITIF D'EGALISATION 
ET DE DECODAGE ITERATIF 
POUR COMMUNICATIONS HAUT-DEBIT 
SUR CANAUX A ANTENNES MULTIPLES 
5 EN EMISSION ET EN RECEPTION 

La presents invention est relative au domaine des communications 
numeriques. . 

Elle conceme plus partlculierement un proc6de et un dispositif 
10 d'egalisatlon et de d6codage Iteratif pour communications haut-d6bit sur 
canaux a antennas multiples en emission et plus particulierement sur 
canaux selectlfs en frequence ^ antennes multiples en emission et 
reception. 

15 DOMAINE TECHNIQUE GENERAL 

Avec le developpement des technologies antennaires, les modeler 
de communication reposant sur le TDMA, le CDMA, POFDM et leurs 
6ventuelles combinaisons sont systematiquement revisit6s et etendus au 
cas de canaux MIMO. Deux grandes classes de codage espace-temps sont 

20 envisagees, qui different dans leurs philosopliies : 

Les transmissions en couches BLAST d6veloppees par les Bell 
Labs ont pour but d'utiliser le muttiplexage spatial pour augmenter le debit 
transmissible sur le canal. II existe deux manieres de coupler codage 
con^ecteur d'erreurs de canal et BLAST, conduisant a deux modeles de 

25 codage espace-temps en couches distlnctes. Dans le premier modele (type 
ST BICM), un codage unique liant les differentes couches est applique sur 
les donnees avant entrelacement spatio-temporel et demultiplexage spatial 
[1, LST-I, p. 1350]. Dans le second modele (type MAC), les donnees sont 
demultiplexees spatialement puis codees en entrelacees independamment 

30 sur chaque couche [1 ,LST-II, p. 1350] ; 

Les schemas de diversite utilisent la diversite spatiale pour mieux 
combattre les effets des evanouissements du canal, et ameliorer les 
performances. On distingue les STTC, qui amenent un gain de codage, et 
les STBC.(et leurs generalisations OD, LDC.:.), qui n'amenerit pas de gain 

35 de 5 Codage et qui sont souvent couples avec les meilleures modulations 
codees en treiilis connues pour le canal gaussien. 
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De r§centes contributions theoriques prouvent que les gains de 
multiplexage spatial et de diversite spatiale sent lies au travers d'une 
relation de compromis. Par exemple, les deux approches pr^cit^es ne 
permettent pas d'augmenter la capacite dans les m§mes proportions pour 

5 un nombre arbitrairement eleve d'antennes d'emission M et d'antennes de 
reception N. On peut verifier que les schemas de diversite a remission ne 
sont optlmaux du point de vue de la theorie de I'information que dans la 
configuration MISO {M, N} = {2,1}, ce qui remet serieusement en cause leur 
pertinence pour les communications haut-debit sur des canaux MIMO a 

10 grand nombre d'antennes. 

Seule I'approche ST-BICM garantif une augmentation lineaire de la 
capacite avec min {M, N}. 

C'est la raison pour laquelle il convlent de privilSgier ce type 
d'emission et de concentrer les efforts sur la conception d'un r§cepteur 
15 adapte garantissant une perfomnance aussi proche que possible des limites 
fortdamentales. A la difference des STBC, les ST-BICM n'assurent pas 
rorthogonaiit6 entre les flux de donnees emis par les diff6rentes antennes. 

Dans le contexte d'une transmission sur canal MIMO selectif en 
frequence, cette strategie de codage appelle done, en reception, une 
20 detection conjointe des donnees pour supprimer a la fols I'interference 
spatiale entre antennes (MAI) et I'interference temporelle entre symboles 
(ISI). 

La fonction de detection conjointe des donnees constitue le point 
delicat et complexe du r§cepteur, et ce d'autant plus que les ST-BICM 
25 necessitent, par essence, un traitement iteratif entre detection et decodage 
de canal pour capturer toute la diversite spatiale. 

PRESENTATION GENERALE DE L'ART ANTERIEUR 
D6coda ae iteratif sur canaux MIMO non selectifs en feouence 
30 Le decodage it6ratif des ST-BICM sur canaux MIMO non selectifs 

en frequence constitue un cas partlculier du probldme d resoudre sous la 
double hypoth^se d'une elimination parfaite de I'ISI et d'une recuperation de 
toute I'energle disponible sur chaque symbole. 

Le prpbleme se reformule alors en ces termes : « suppression de la 
35 MAI sur un canal MIMO plat en presence de baiit colore spatialement ». II 
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necessite le calcul des APP sur les bits des symboles des vecteurs transmis 
a chaque Instant (statistique sur les donnees a transmettre au decodeur). 
Deux methodes ont ete retenues dans la Iltt6rature : 

Calcul exact des APP base sur une liste exhaustive (MAP). La 
5 complexite de cette methode est prohibitive, puisque exponentielle en le 
nombre d'antennes emettrices M et polynomlale en Tordre Q des 
constellations employees. 

Calcul approche des APP a partir d'une llste r6duite de 
vecteurs candidats gen§ree par un algorlthme de decodage par sphere [2]. 
10 L'algorithme de d6codage par sphere pr6sente au moins deux avantages : il 
preserve le crit^re MAP et possede une complexite en O (M3), insensible a 
i'ordre Q des constellations employees. On precise que Talgorithme de 
decodage par sphere peut §tre vu comme une variante de l'algorithme de 
decodage sequentiel de Fano. 
15 - Calcul approche des APP a partir d'une liste reduite de 

vecteurs candidats gener6e par un algorlthme sequentiel employant une 
pile (stack algorithm). 

L'algorithme de decodage par sphere a ete intensivement etudie 
pour les canaux MIMO non selectifs en frequence. Sa generalisation au cas 
20 des canaux MIMO selectifs en frequence est non triviale. 

Une approche brutale serait de rappliquer par blocs. Cette approche 
est peu interessante car elle augmente artificiellement la dimension de 
I'espace de recherche (et done la complexite) et introduit un decodage 
redondant des symboles. 
25 L'algorithme decrit dans [2] est derive pour des canaux MIMO non 

selectifs en frequence. Sa generalisation au cas des canaux MIMO selectifs 
en frequence est un probl^me ouvert. 

La performance de ces deux algorithmes est essentiellement 
detenminee par la qualite de la liste de candidats generee, 
30 pr^ferentiellement de taille fixe. II est primordial qu'elle comporte le meilleur 
candidal decodant un bit 1 et le meilleur candidat decodant un bit 0 a toute 
position. 

Decodaoe iteratrf sur canaux MIMO selectifs en freouence 
Plusieurs types de detecteurs inseres dans des structures iteratives 
35 ont ete proposes. 
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Detecteur a entrees et sorties ponderees reallsant le critere du 
MAP base sur ralgorithme BCJR. La complexite de ce recepteur, en 
O(QLxM), est ciairement prohibitive pour des modulations d'ordre eleve et 
des canaux MiMO a grand nombre d'entrees et memoire elevee. 

5 - Detecteurs a entrees et sorties ponderees sous-optimaux 

bases sur des algorithmes de recherche efficaces dans des treillis fortement 
reduits en nombre d'etats [3], et, le cas echeant, precedes de filtres a 
minimum de phase. Cette approche est encore limitee par sa complexite, en 
O(QLrxM) ou Lr represente la longueur de contrainte reduite. En particulier, 

10 Temploi de constellations a grand nombre d'etats est exclu. 

Les detecteurs non lineaires construits a partir d'une annulation 
lineaire iterative de I'interference ([4] dans un contexte de detection multi 
utilisateur). La transposition de ces d§tecteurs dans un contexte MIMO est 
en generate bases sur Tanalogie entre antennes et utilisateurs diff^rents. 

15 Leurs implementations n6cessitent done un filtre par antenne et leurs 
performances sont limitees par I'approche disjointe MMSE par antennes. 
Cette approche presente I'avantage d'avoir una complexite polynomiale en 
nombre d'antennes d'emission. 

20 PRESENTATION DE L'lNVENTION 

Un but de Tinvention est de proposer un recepteur avance pour les 
transmissions numeriques haut-debit pour canaux selectifs en frequence a 
antennes multiples en emission et reception, qui soit peu complexe - et 
notamment qui ne necessite pas des puissances de calcul importantes- 

25 tout en permettant de traiter Tinterference spatiale entre antennes et 
interference temporelle entre symboles. 

A cet effet, I'invention propose un dispositif d'egallsation et de 
decodage iteratif pour communications haut-debit sur canaux selectifs en 
frequence a antennes multiples en 6mission et en reception, comportant un 

30 egaliseur a retour de decision qui regoit les donnees des differentes 
antennes de reception, cet egaliseur comportant un filtre avant, ainsi qu'un 
filtre arriere recursif aliments par une reconstitution ponderee des donnees 
calcuiee a partir de la sortie d*un decodeur lui-m§me alimente par des 
moyens de prise de decision, le dispositif comportant en outre des moyens 

35 pour soustraire la sortie dudit filtre arriere aux donnees en sortie du filtre 
avant, les donnees ainsi soustra'ites etant transmises en entree des moyens 
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de prisfe de decision avec la sortie du decodeur, les moyens de prise de 
decision generant en outre une statistique transmise a destination d'un 
decodeur de canal a entrees et sorties ponderees, lesdits moyens de prise 
de decision prenant en compte la correlation spatiale du bruit en sortie des 
5 moyens de soustraction, lesdits moyens de prise de decision et le decodeur 
etant separes par des precedes d'entrelacement spatio-temporel au niveau 
binaire, 

caract^rise en ce que le filtre avant et arrlere sont adapte de fagon 
iterative pour minimiser Terreur quadratique moyenne en sortie du 
10 soustracteur. 

Un tel dispositif est avantageusement complete par les differentes 
caracteristiques prises seules ou selon toutes leurs combinaisons 
techniquement possibles : 

les moyens de d6cjsion en sortie des moyens de soustraction 
15 de r6galiseur sont de type blanchiment spatial suivi d'un decodage par 
sphere ; 

les moyens de decision en sortie des moyens de soustraction 
de Tegaliseur sont de type SIC/PIC (suppression d'interference spatiale 
residuelle en sortie des moyens de soustraction de Tegaliseur, en serie 
20 et/ou en parallele) ; 

le blanchiment spatial est mis en oeuvre en sortie des moyens 
de soustraction de Tegaliseur ; 

le blanchiment spatial est mis en oeuvre par les moyens de 
prise de decision ; 

25 . - le blanchiment spatial est mis en oeuvre par les filtres avant et 

arriere ; 

a partir d'une certaine iteration, le filtre avant est un filtre 
adapte, le filtre arriere le meme filtre adapts moins le coefficient central. 

L'invention decoule d'une approche originate radicalement differente 

30 de celles considerant les antennes comme des utilisateurs differents. En 
effet, le signal transmis est vu comme une modulation T-dimensionnelle, T 
etant le nombre d'antennes en emission. Ce changement de point de vue a 
des consequences considerables sur le design du recepteur. Le recepteur 
se decompose en un egaliseur lineaire iteratif, considerant k son entree une 

35 modulation vectorielle T-dimensionnelle convoluee a un canal selectif en 
frequence, suivi d'un detecteur de modulation T-dimentionnelle (par 
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example sphere d^codeur) capable de generer des informations souples 
pour le decodeur de canal. Cette approche presente deux avantages par 
rapport a Tart anterieur. Elle est moins complexe car elle necessite un 
unique filtre (vectorielle) et est plus perfomnante car elle permet le choix du 

5 critere de la detection vectorielle. 

Uinvention conceme egalement un systeme de communication 
haut-debit sur canaux selectifs en frequence a antennes multiples en 
emission et en reception, caracteris6 en ce qu'il comporte en reception un 
dispositif d'egalisatlon et de decodage du type precrte. 

10 Avantageusement notamment, les moyens d*emission sont de type 

ST-BICM. 

En outre, Tinvention propose un proc6de d'egalisatlon et de 
decodage iteratif pour communications haut-debit sur canaux s§lectifs en 
frequence ^ antennes multiples en emission et en reception, mettant en 

15 oeuvre un egaliseur d retour de decision qui regoit les donnees des 
diff6rentes antennes de reception, cet egaliseur comportant un filtre avant, 
ainsi qu'un filtre arriere recursif aliments par une reconstitution ponderee 
des donn§es calcul6e a partir de la sortie d'un decodeur lui-meme alimente 
par des moyens de prise de decision, le dispositif comportant en outre des 

20 moyens pour soustraire la sortie dudit filtre arriere aux donnees en sortie du 
filtre avant, les donnees ainsi soustraites etant transmises en entree des 
moyens de prise de decision avec la sortie du decodeur, les moyens de 
prise de decision generant en outre une statistique transmise a destination 
d'un decodeur de canal a entrees et sorties ponderees, lesdits moyens de 

25 prise de decision prenant en compte la correlation spatiale du bruit en sortie 
des moyens de soustraction, lesdits moyens de prise de decision et le 
decodeur etant separes par des precedes d'entrelacement spatio-temporel 
au niveau binaire, 

caracterise en ce que le filtre avant et anri§re sont adaptes de fagon 

30 Iterative pour minimiser I'erreur quadratique moyenne en sortie du 
soustracteur. 

PRESENTATION DES FIGURES 

D'autres caracteristiques et avantages de i'invention ressortiront 
35 encore de la description qui suit, laquelle est purement illustrative et non 
limitative et doit §tre lue en regard des dessins annexes, sur lesquels : 
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La Figure 1 illustre le concept de la technique V/BLAST ; 
- La Figure 2 illustre le modele g6n6ral de communication 
connu sous le nom de ST-BICM ; 

La Figure 3 illustre Tarchitecture d'un recepteur iteratif 
5 conforme a un mode de realisation de I'invention ; 

La Figure 4 Illustre I'architecture d'un recepteur iteratif qui peut 
§tre utilise a partir par exemple de la deuxi§me iteration ou des iterations 
suivantes. 

10 DESCRIPTION D'UN OU PLUSIEURS MODES DE REALISATION 

DE L'INVENTION 

Rappels sur la technique BLAST et le codaae ST BfCM 

On a represents sur la figure 1 le schema general d'une architecture 

BLAST. 

15 Les donndes codSes dans un dispositif d'encodage vectoriel EV 

unique, sont transmises, apres entrelacement spatio-temporel et 
d6multlplexage spatial, sur une pluralite d'antennes d'emission TX ; a I'autre 
extremite des canaux MIMO, les signaux emis sont regus par une plurality 
d'antennes de reception RX qui les renvoient sur des moyens de dScodage 

20 DE (canaux MIMO) en sortie desquels on recupSre les donnees. 

Dans ce qui est d§crit ci-apres, on se place dans le cas de canaux 
MIMO selectifs en frequence. La presence d'interferences entre symboles 
augmente la complexite du recepteur. 

La figure 2 illustre quant a elle le module general de communication 

25 connu sous le nom de ST BICM. 

A remission, les donnees sont traitees par un code correcteur de 
canal 1 (codes convolutifs, turbo-codes, codes LDPC...), puis font I'objet 
d'un entrelacement appliqu§ au niveau binaire 2, d'un demultiplexage 
spatial 3, et, pour chaque couche, d'un precede de modulation 4-1 ^ 4-M. 

30 Apres un filtrage de mise en forme 5A a 5M, ies donnees ainsi moduiees 
sont transmises a des antennes emettrices 6-1 d 6-M. 

Description des traltements d'emission et de reception 
A remission, les donnees font I'objet d'un traltement ST-BICM. 
35 A cet effet, on met en oeuvre les etapes suivantes : 

Reception d'un flux de donnees numeriques a un debit donne ; 
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Application d'un codage correcteur de canal 1 generant un flux 
de donnees numeriques codees ; 

Entrelacement des donnees numeriques codees par un 
entrelaceur 2 operant au niveau bit ; 
5 - Demultiplexage spatial 3 du flux de donnees numeriques 

codees entrelacees pour creer une pluralite de flux distlncts de donnees 
numeriques codees appeles couches emettrices (wflux distincts ou couches 
emettrices) ; 

Puis, pour chacun des M flux distincts de donnees numeriques 
10 codees : 

Modulation du flux de donn6es numeriques codees 
entrelacees 15 conformement d un sch§ma de modulation (4-1 ^ 4-M) pour 
obtenir un flux de symboies modules ; 

Filtrage du flux module (5-1 a 5-M) ; 
15 - Transmission du flux module sur une antenne d'emlssion 

propre (6-1 d 6-M). 

Le recepteur du systeme de transmission est illustre sur la Figure 3. 
II comprend un egaliseur d retour de decisions defini par un filtre 
avant 9, un soustracteur 10, un algorithme prise de decisions 11 et un filtre 
20 arridre recursif 1 2. 

Le filtre arriere recursif 12 est alimente par une reconstitution 
ponderee des donnees calculee a partir de la sortie d'un decodeur13 lul- 
meme alimente par Talgorithme de prise de decision 1 1 . 

Le filtre avant 9 et le filtre. arriere 12 sont detemnines de faQon 
25 iterative pour minimiser I'erreur quadratique moyenne (MSE) de Tegaliseur 
MIMO, c'est-a-dire pour minimiser Terreur en sortie du soustracteur 10. 
lis utillsent a cet effet une estimation Initiale 8 du canal MIMO. 
Le traftement mis en oeuvre est par exemple du type de celui 
detaille en annexe II. 
30 On notera que, dans le cadre ce traltement. on travaille sur une 

estimation vectorielle de I'erreur residuelle et du bruit, les filtres avant et 
arrldre etant calculus par blocs pour minimiser ce vecteur. 
Ce traltement vectoriel permet un calcul simplifie. 
Les moyens de prise de d6clsion11 regoivent en entree les 
35 donnees en sortie du soustracteur 10, ainsi que la sortie du decodeur 13. 
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L'algorithme qu'ils mettent en oeuvre peut etre de different type. II 
peut notamment etre de type SIC/PIC (suppression d'Interference spatlale 
residuelle en sortie des moyens de soustraction de Tegaliseur, en serie 
et/ou en paralleie) (Annexe III). 
5 En variante. il peut etre un algorithme de decodage par sphere. 

Un tel algorithme possede une complexite en 0(M^) - oCi M est le 
nombre d'antennes d'emission - ind^pendante de Q. 

Ceci rend possible I'utilisation de modulations a grand nombre 
d'etats en vue d'augmenter le d§bit. 
10 La sortie des moyens 11 de prise de decision est transmise a un 

proc6d6 de desentrelacement spatlo-temporel 14. qui est mis en oeuvre au 
niveau binaire et qui est interpose entre lesdits moyens de prise de 
decision 1 1 et le decodeur 13. 

La sortie du decodeur 1 3 est une probabilrte sur bits. 
15 Cette probability est transmise d'une part en entree de 

l'algorithme 1 1 de prise de decision et d'autre part, apr§s entrelacement 
spatio-temporel (proc6d6 15). S un traitement 16 de reconstltution ponder^e 
des donnees. 

Les donnees ponder§es ainsi reconstituees sont celles qui sont 
20 transmises en entree du filtre aniere 12. 

La sortie du filtre amere 12 conBspond a une r^gyn^ratioh 
ponderee des donnees. 

Par ailleurs, on notera que I'en-eur con-espondant a Tinterf^rence 
residuelle et au bmit qui est injects en entree du soustracteur 10. est 
25 coloree ^ la fois temporellement et spatialement. 

SI la con-elation temporelle a peu d'impact sur les traitements de 
regaliseur, la correlation spatlale joue un role essentiel. 

C'est pourquoi, dans un mode de realisation, on realise un 
blanchlssement spatial par factorisation de Cholesky. 
30 Ce blanchiment spatial est avantageusement mis en oeuvre en 

sortie du soustracteur 10. 

Dans une variante en partlculler, il peut §tre mis en oeuvre par 
l'algorithme 1 1 de d§codage par sphere lui-meme. 

II peut egalement §tre mis en oeuvre en 6tant Integra aux filtres 
35 avant et aniere 9 et 1 2. 
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Enfm, dans une variants de mise en oeuvre possible, le filtre avant 9 
est avantageusement remplace par un filtre adapte a partir d'un certain 
nombre d'iterations donnees, par exemple a partir de la seconde iteration. 

Comme le filtre arriere se deduit directement du filtre avant, il 
5 resulte done, dans ce cas, de la convolution du filtre adapte avec le canal 
moins le coefficient central. 

C'est ce qu'illustre la figure 4, sur laquelle le filtre avant 9 est un 
filtre adapte MF, tandis que le filtre anidre 12, denomme Bmf, est ie filtre 
adapts convolue avec le canal moins le coefficient central. 
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ANNEXE I 

LISTE DES ACRONYMES 

TDMA : Acces multiple a repartition dans ie temps (Time Division 
5 Multiple Access) 5 CDMA : Acces multiple par repartition dans Ie code 
(Code Division Multiple Access) 

OFDM : Multiplexage par repartition orthogonale de la frequence 
(Orthogonal Frequency Division Multiplex) 

MIMO : Entree multiple sortie multiple (Multiple Input Multiple 

10 Output) 

SISO : Simple entr6e simple sortie (Single Input Single Output) 
BLAST : Couche spatio-temporelle des laboratoires Bell (Bell Labs 
. Layered Space Time) 

ST-BICM : Modulation codee spatio-temporelle avec entrelacement 

15 binaire Space-Time Bit-Interleaved Coded Modulation) 

STTC : Code trelllls spatlo-temporel (Space-Time Trellis Codes) 
STBC : Code bloc spatio-temporel (Space-Time Block Codes) 
MAI : Brouillage acces multiple (Multiple Access Interference) 
ISI : Interference inter symbole (Inter Symbol Interference) 

20 APP : Probabilite A Posteriori (A Posteriori Probability) 

BCJR : Bahl, Cocke, Jelinek. Raviv (algonthme de) (Bahl. Cocke. 
Jelinek, Raviv (algorithm of)) 

• MF : Filtre adapte (Matched Filter) 
WMF : Filtre adapte blanc (Whitening Matched Filter) 

25 MAP : Maximum A Posteriori (Maximum A Posteriori) 

MMSE : Erreur quadratique moyenne minimum (Minimum Mean 
Square Error) 

IC : Annulation d'interference (Interference Cancellation) 
Soft-IC : Annulation d'interference souple (Soft-Interference 
30 Cancellation) 

SIC: Annulation d'interference serie (Serial Interference 

Cancellation) 

PIC : Annulation d'interference paralldle (Parallel Interference 
Cancellation) 

35 LDPC : Controle de parite de faible density (Low Density Parity 

Check) 
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ANNEXE II 

MODELE DE COMMUNICATION 

On considere un canal P-bloc a entrees et sorties multiple a 
5 T antennes d'emission et R antennes de reception selectif en frequence de 
memoire M. 

A. Modulation codee d'entrelacement binaire spatio-temporel 
Soit C un code lineaire de longueur Nc et de rendement pc sur F2 admettant 
10 un vecteur de donnee u e ^2^^^ et prdduisant un mot code c ^ Fz^. Nous 
supposons que le rendement pc indue des bits de queue si des codes 
conventionnels sont employes. Le mot cod6 entre dans un entrelaceur 
binaire bien congu n, dont la matrice de sortie A e F2^^*^*" est segmente en 

P matrices A*' e F2^^*^^ . P = 1 P- Les colonnes de la matrice A*' sont 

15 des vecteurs a^ [n] e Fz^, n = 1, .... L, que Ton nomme « vecteur 
d'etiquette de symbole », contenant T sous vecteurs a^[n] e F2^, t = 1, 
T (un par canal d'entree), avec les composant binaires empll6s a^<tj>[n], ... 
a^<tj>[n], ou <tj> repr§sente I'index (t - 1) q+j. Dans chaque matrice A^, 
tous les vecteurs a^[r\] sont modules par un modulateur sans memoire D- 
20 dimensionnel sur un ensemble de signal A cz c"^ de cardinality |A| = 2^^ a 
travers le codage bit a symbole p : {0, 1}^^ A. Sans perte de generalite, 
nous supposons que A est un produit simple de T constellations complexes 
C de meme ordre 2^ et que par constellation, un codage de Gray est utilise. 
Apres transformation du signal, ce processus de modulation codee peut 
25 etre aussi considere comme un schema de modulation codee spatio- 
temporelle ou tous les symboles complexes de mot code RxLxT sont 
divises en P matrices Xp e C^'^S P = 1, P, dont les colonnes x^ [n] e C^, 
n = 1, L sont nommees « vecteur de symbole de constellation ». Le 
« vecteur d'etiquette de symbole » a" [n] peut etre obtenu de x^ [n] par une 
30 simple inversion de codage bit a symbole En tombant dans la classe 
generale des codes spatio-temporels, on appelle cette architecture 
Modulation Codee Spatio-Temporelle avec Entrelacement Binaire 
(STBICM) [3]. D'apres Thypoth^se de filtrage ideal de Nyquist en bande 
Ilmit6e, le rendement de transmission (en bits par canal utilise) est : 

35 

Un schema de principe relatif a ce modele de communication est montre 
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figure 1 . 

B. Canal MIMO 

5 Soit HP e c^''^"^^''^^ le canal MIMO pour le bloc p, oD p = 1, P et H = 
{H^} Tensemble de tous les canaux. Soient egalement e c'^'*^ et e 
C'^*'^ respectivement le « matrice de symbole de constellation » et la 
« matrice de sortie de canal ». La sortie vectorielle de canal en bande de 
base equivalente a temps discret [n] e au temps n = 1, .... L peut etre 

10 ecrite comma : 

[n] ^ (H ( i }] X? H + N * I] W X? (n - *1 + [n] (2) 

Ou [n] ^ C"^ sent le « vecteur de symbole de constellation » transmis aux 
temps n, I'energie de chaque composante etant egale a ax^ H*^ [k] ^ C^""^ 
est la matrice du nombre de taps k de la reponse impulsionnelle du canal, 

15 [n] e est le vecteur de bruit complexe additff. Les vecteurs de bruit 
complexe addltif [n] sont supposes de moyenne nulle, independant 
Identiquement dlstribu6s, et Gaussien complexes a symetrie circulaire, et 
suivent done le pdf N (0, o^l). Le canal H", constant pendant la duree du 
bloc correspondant, a une reponse impulsionnelle finie (FIR) de longueur 

20 M+1, dont les taps de symbole d'espace [0] , ...,H^ [M] sont les matrices 
aleatoires complexes RxT, de moyenne nulle et dont la puissance moyenne 
satisfait les contraintes de normalisation : 

B Jdiag |]CH? {k] IF (Af I j « n (3) 

dans le cas d'un systeme de puissance 6gale. L'operateur « t » correspond 
25 a Toperateur de transposee conjugu§e. La memoire de canal egale pour 
tous les liens RxT possibles est une hypothese raisonnable 6tant donne 
que le nombre de composant individual multivoies est dicte de maniere 
pr§dominante par de larges structures et des objets reflechissants. 
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Bloc iteratif d'ahnulation d'interf6rence MIMQ 
A. Principe et ordonnancement 

5 Le bloc iteratif M\MO annulateur d'interference traite cliaque bloc de donn6e 
regus Y^, p = 1, .... P , sepai^ment, en utilisant I'lnfonmation probabillste 
al^atoire renvoy^e par le decodeur de sortie . Tous les signaux et 6l6nrients 
d^crits sont done Indexes par bloc. Nous supprlmons quelque fols cette 
dependence pour simplifier les notations. Durant chaque iteration I, un filtre 

10 avant lineaire F\ applique a cliaque symbole regu de vecteur y [n], produit 
le signal de vecteur y* [n]. Puis une estimation construite de mani^re 
appropri§ e' [n] du MAI et du ISI alterant x [n] est soustraite du signal de 
vecteur y" [n] pour produire z" [n] : 

15 L'estimation du vecteur e' [n] provient de la sortie du filtre arriere B' excite 
par le vecteur x* [n] de tentative de decision souple sur le vecteur x[n], etant 
donne la connaissance prealable disponible (distribution de probabilite 
extrinsdque) ^ I'lt^ration I - 1 . 

Nous amvons maintenant au coeur du turbo-egaliseur MIMO, la derivation 
20 des flltres a r^ponse impulsionnelle finie MIMO ayant et arriere. 

B. Calcul des filtres avant et arrifere 

Nous commengons par quelques manipulations de base sur la sortie 
25 instantan6e y'[n] du filtre avant F* : 

y'W^ Z P'WylM-iJ (22) 

Puisque : 

M 

ytrtl = 5^Hlfclx[rt-A!]+w[nl (23) 

fc=0 

Chaque 6chantlllon y'[n] peut §tre d6velopp§ comme suit : 
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<25) 



i|!4 Af 

y*l"H J2 5;F'KHtft]x[n-»-fe]+ 52 F'lifwt^i-i] (24) 

I!>eO isi—hf^ 

Une notation matricielle Squivalente pour la convolution de deux filtre He et 
F' est: 

5 OD: 

est le filtre avant d'ordre Lf = Lfi + Lf2+1 . ou : 
H[0] H[i] < - HlJkf] 

H(0] H[ll • • H[M] 

H[0J H[l] H[JI<| 



(26) 



10 



15 



H[0] H[ll HJMI 

est la matrice de canal de bande-diagonale de Toeplitz et ou 

x[n + Lj^\ 



3teN"[ x[nl ec^t^^*^ 

est le vecteur transmis. 

En introduisant le filtre combine : 

On obtient I'expresslon finale : 

Avec Lgi = Lpi, Lg2 = Lf2 + M et Lg = Lf + M. 



<2S) 



(30) 



(31) 



Par analogie, on definit le filtre arriere comme suit 
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[b'[-Lb,\ ■■' 0 -•-B'[iftl] eC'^- (32) 
La sortie du correcteur de brouillage peut etre 6crite : 

2' [n] = ^ G' m ^[n-k\-J2^'[k]^[n-ki-\-Y,¥'[k]w\^- k] (33) 
fe fc fc 

Le vecteur de bruit et de brouillage r6siduel \/[n] (d6nomm6 aussi vecteur 
5 d'erreur de I'egaliseur Ml MO) est exprim§ comme suit : 

V[n] = as* [rt] - G' lOjx^il 

- T.^' M*t« - - ebuj^s^ [rt - &] + x;^' Sfe] w [« - ^^l 

= 5;;fi'[ifc]x{n-A^-XlB'[*:lS'b»-*l + I2^^W«'tn-&l (34) 



OD : 

G' - [G' [-La,] - . 0 - . G* {LoJi\ (35) 

10 Les flltres F' et B' sent calcules par blocs pour minimlser I'erreur quadratlque 
moyenne de I'egaliseur MIMO, sous la contrainte G'[0] = I, un probleme qui, 
dans le cas MIMO. peut etre formula de manidre compact oomme suit : 

fF^B^I s= arg inin trE-fV N v^fn - fef ) 

=s arg inin trfKl} 

La minimisation MSE est executee en deux 6tapes successives et donne 
15 finalement : 

B'=fif (37) 
P' = (EAHj*'"'HeEL)"*EAHl*'~' (38) 

OCi: 



Ea = 1- 0 — I-.. 0- -l€C=^^ (39) 



et ou : 
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= K*^ - H^rHj + «r»l] (40) 

avec : 

n=l 

Nous insistons sur le fait que le baiit plus le brouillage restduel a la sortie du 
.5 soustracteur est correle spatialement et temporellement. La correlation 
temporelle n'a pas d'impact reel sur les traitements suivant, la correlation 
spatiale joue un r6le ici essentiel. Le bruit blanc spatial peut etre aisement 
obtenu par une simple factorisation de Cholesky de la matrice de correlation 

K'j». 

10 Proposition 2 : Puisque la matrice de correlation k'j, est definie positive, la 
factorisation de Cholesky est toujours applicable. Sachant cela : 
K^«LL» (^) 

Ou L est une matrice triangulaire inferieure. Ensuite, en appliquant f'' = L'** 
F' et B'' = L'"* B' comme fiitre avant et filtre arriere (au lieu de f' et B'), la 
15 matrice de con-elation de brouillage plus bruit egale ^ la matrice identite. 
Preuve : La preuve est immediate. 

= L-^ [K - "4) ^JS'* -f o^F"*] - ir^K = I <43) 

Vailde Talgorithme propose par la simulation de la figure 3. 
Dans la suite, nous considerons toujours la sortie de I'egaliseur avec 
20 blanchlment spatial de sorte que : • - 

z'[fi]:=G'[0]x[n] + VM 

Avec Kv' = I et G'[01 = L-^ 

II «st int§ressant de remarquer que le critere MMSE correspond a une 
maximisation du SNR tr{G'[0]*Kv'G'[0]}. Le SNR maximal donne pour le filtre 
25 adapte est atteint pour = cxx^- 
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C. D6codeur APP exhaustif 

Nous definissons ici toutes les quantit^s probabilistes Schang^es au cours 
5 des iterations. La partie detection de I'egaliseur MiMO a pour but de foumir 
les informations extrinsdques sur les bits des symboles de la matrice A 
selon le nouveau moddle de communication (10) ou G' [0] joue le role d'un 
canal plat MIMO equivalent T x j et v* [n] est I'interfdr^nce r^siduelle plus 
bruit blanchi. Les rapports logarithmique APR sur tous les bits des symboles 
10 a<tj>[n] sont definis comme suit : 



15 



(H) 



Par simple marginalisation, on obtient : 



OU : 



Les rapports logarithnfiiques APP peuvent dtre developp§s comme : 



20 En supposant un entrelacement spatio-temporel parfait autorise de diffuser 
les probabilites pr6c6demment Jointes dans le produit de leur composant a 
chiifre marginal : 

En supposant que [n] a ete blanch! spatialement par un moyen 
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quelconque (exempie factorisation de Cholesky), la metrique Euclidienne 
peut §tre utilises dans revaluation de la probability : 

p{^'[4\^) cxcxp|-^||z'H-GMojx||"} (16) 

Dans le cas general, la matrlce de correlation K'v du vecteur de bruit doit 
5 etre prise en compte dans revaluation de la probabilite : 

j»(z'[»]|x)ae5q>{-(z*[n]-G'IO]x)*KirM2'N-G'LOlx)} (17) 

Comme il est etabii dans ie domaine des turbo-detecteurs, les taux de 
probabilite extrinseque logarithmique sont deduit selon la formule : 

10 Tous les echantlllons d'informatlon extrinseque logarithmique du detecteur 
MIMO sont collectes a partir de tous les blocs p = 1, .... P et rearranges, 
apres le desentrelacement spatio-temporel n'\ en un vecteur d'observation 
simple ^ ^ e R^, sur la base duquel le decodeur de sortie foumi les rapports 
de probabilite extrinseque logarithmique sur tous les bitis de mots code : 



Pr 



15 



20 



Apres le desentrelacement spatio-temporel le vecteur est dif^jse dans 
P matrices n'"" de teux de probabilite precedente. logarithmique, un pour 

chaque bloc de donnees A**. Pour chaque bloc p = 1 P, les entrees de 

sont donnees par : 

Ceci permet de reecrire le vecteur x" [n] de prise de decision soft comme : 

^N=iE-n(i+(^s>(-)-0*^(^)} (21) 
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D. Detection MIMQ Quasi-ODtimal oar i'interm^diaire du decodeur a sphere ^ 
lists modifiee 

La detection iVIIMO optimal ne peut pas etre retenue pour un sc^narii de 
5 communication haut d6bit. En effet, {'^numeration de point des sous- 
ensembles entiers de constellation n*°^<tj> et n*^*<tj>. dont la cardinalite 
varie dans 0(2''^), peut rapidement devenir accablante pour des 
modulations d'ordre sup§rleur et/ou un grand nombre d'antennes de 
transmission. Una analyse soigneuse des valeurs de probability indique 
10 qu'un grand nombre d'entre elles sont negligeables. Par consequent, 
comme contribution significative de cet article, nous suggerons de 
remplacer renum6ration approfondie de point par une enumeration de point 
des sous-ensemble L^°*<t.j> et L*'''*<tj> de taille fortement reduite, sous 
ensemble aussi appel6s liste, qui contiennent uniquement las probabilites 
15 non n6gligeables. Les rapports de probabilite extrinseque logarithmique 
deviennent ators : 

En temies g^om^triques, contiennent des points du treillis dans une sphere 
20 de rayon r centr6e sur un point bien choisi (par exemple le point ML sans 
contrainte ou le point ML lui-m§me). Des versions modifi6es de decodeur a 
sphere ont 6te employees avec un certain succ§s pour exhiber ces listes 
[2]. Le choix du rayon de sphdre r pilots la quality et la complexite du 
. detecteur MIMO correspondant. 

25 

E. Alaorithme de detection PIC 

Le decodeur a sphere peut egalement etre considere comme trop 
complexe. Una mani^re possible d'alleger encore la complexite du 
30 detecteur est de generer I'APP log independamment par dimension. La 
variable de decision pour le composant x est, a i'iteration I : 
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OD Qk est la coionne k de la matrice Gl[0]. Les rapports APP log ^'<tj>[n], j= 

1 q sont alors calcules en supposant que Ht[n] est une variable 

gaussienne de moyenne Xt, et de variance estim6e : 

5 

Le premier detecteur d'iteration peut etre ameliore en utllisant une 
technique SIC de la fagon suivante : 
• Initialisation 

Classifier les composants par ordre decroissant en fonction de leur 
10 rapport signal sur parasite (SIR), par exemple : 



avec 



SIR(xt) > SIR(X2) >"> Sm{xT) 



Recursivite 
15 Pour t =1 a T 
Calcul de : 



Calcul du rapport APP log ^<g>[n] , j= 1 q en supposant que r^Jn] est une 
variable gaussienne de moyenne Xt, et de variance estimee : 



1 

L 

20 

Calcul I'estimation MMSE de connaissant Ht[n] 
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Xt[n]^E{xtln]\rl[n]} 
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REVENDICATIONS 

1. Dispositif d'^galisation et de decodage iteratif pour 
communications haut-debit sur canaux selectifs en frequence a antennes 

5 multiples en emission et en reception, comportant un egaliseur a retour de 
decision qui regoit les donnees des differentes antennes de reception, cet 
egaliseur comportant un filtre avant (9), ainsi qu'un filtre anriere (12) recursif 
alimente par une reconstitution ponderee des donnees caiculee a partir de 
la sortie d'un decodeur (13) lui-meme alimente par des moyens de prise de 

10 decision (11), le dispositif comportant en outre des moyens pour souslraire 

(10) la sortie dudit filtre anidre (12) aux donn6es en sortie du filtre avant (9), 
les donn6es ainsi soustraites etant transmises en entree des moyens de 
prise de decision (11) avec la sortie du d6codeur (13). les moyens de prise 
de decision (11) g6n6rant en outre une statistique transmise d destination 

15 d'un decodeur de canal § entrees et sorties pond6r6es, lesdits moyens de 
prise de decision (11) prenant en compte la correlation spatiale du bruit en 
sortie des moyens de soustraction (10), lesdits moyens de prise de decision 

(11) et le d§codeur (13) 6tant sSparSs par des proc6des d'entrelacement 
spatlo-temporel au niveau binaire, 

20 caract6rise en ce que les filtres avant (9) et arriere (12) sont 

adaptes de fagon iterative pour minlmlser I'enBur quadratique moyenne en 
sortie du soustracteur (1 0). 

2. Dispositif selon la revendication 1 , caracterise en ce que les 
moyens .de decision (11) en sortie des moyens de soustraction (10) de 

25 I'egaliseur sont de type blanchiment spatial suivi d'un decodage par spliere ; 

3. Dispositif selon la revendication 1 , caracterise en ce que les 
moyens de decision (11) en sortie des moyens de soustraction (10) de 
I'egaliseur sont de type a suppression d'interf§rence spatiale r§siduelle en 
sortie des moyens de soustraction (10) de Tegaliseur, en serie et/ou en 

30 parall§le(SIC/PlC). 

4. Dispositif selon I'une des revendications pr6cedentes. 
caracterise en ce que le blanchiment spatial est mis en oeuvre en sortie des 
moyens de soustraction (10) de I'egaliseur. 

5. Dispositif selon la revendication 4. caract6ris6 en ce que le 
35 blanchiment spatial est mis en oeuvre par les moyens de prise de decision 

(11). 
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6. Dispositif selon la revendication 2, caracteris6 en ce que le 
blanchiment spatial est mis en oeuvre par les filtres avant (9) et arriere (12). 

7. Dispositif selon Tune des revendications precedentes, 
caracterise en ce que a partir d'une certaine iteration, le filtre avant (9) est 

5 un filtre adapte. 

8. Systeme de communication haut-debit sur canaux selectifs en 
frequence a antennes multiples en emission et en reception, caracterise en 
ce qu'il comporte en reception, un dispositif d'egalisation et de decodage 
selon Tune des revendications precedentes. 

10 9. Systeme selon la revendication 7, caracterise en ce que les 

moyens d^emission sont de type ST-BICM. 

10, Procede d'egalisation et de decodage it6ratif pour 
communications haut-debit sur canaux selectifs en frequence ^ antennes 
multiples en Emission et en reception, mettant en oeuvre un egallseur a 

15 retour de decision qui regoit les donnees des diff^rentes antennes de 
reception, cet §gali$eur comportant un filtre avant (9), ainsi qu'un filtre 
arriere (12) r§cursif alimente par une reconstitution ponderee des donnees 
calcujee d partir de la sortie d'un d§codeur (13) lui-m§me alimente par des 
moyens de prise de decision (11), le dispositif comportant en outre des 

20 moyens (10) pour soustraire la sortie dudit filtre arriere (12) aux donnees en 
sortie du filtre avant (9), les donnees ainsi soustraites etant transmises en 
entree des moyens de prise de decision (11) avec la sortie du decodeur 
(13), les moyens de prise de decision (11) generant en outre une statistique 
transmise a destination d'un decodeur de canal a entrees et sorties 

25 ponderees, iesdits moyens de prise de decision (11) prenant en compte la 
correlation spatiale du bruit en sortie des moyens de soustraction (10), 
Iesdits moyens de prise de decision (11) et le decodeur (13) etant separes 
par des precedes d'entrelacement spatio-temporel au niveau binaire. 
caracteris6 en ce que le filtre avant (9) et arriere (12) sont adaptes de fagon 

30 iterative pour minimiser Terreur quadratique moyenne en sortie du 
soustracteur(IO). 

11. Procede selon la revendication 10, caracterise en ce que le 
filtrage recursif de i'egalisation a retour de decisions est egalement adapte 
de fagon Iterative pour minimiser Terreur quadratique moyenne en entree de 

35 Fegallseur. 
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.12. Procede selon I'une des revendications 9 ou 10, caracterise 
en ce que Von met en CBUvre un blanchiment spatial sur les donnees 
soustraites. 

13. Procede selon Tune des revendications 9 a 1 1 , caracterise en 
5 ce que ralgorithme de prise de decision (11 ) est de type decodeur par 

sphere. 

14. Procede selon les revendications 11 et 12 prises en 
combinaison, caracterise en ce que le blanchiment spatial est mis en ceuvre 
par ralgorithme de prise de decision (11). 

10 15. Procede selon Tune des revendications precedentes, 

caracterise en ce qu'& pari:ir d'une cerlalne iteration, le filtre avant (9) est un 
filtre adapte. 



wo 2005/025118 



PCT/FR2004/002104 



1/2 




FIG.1 




FIG.2 



wo 2005/025118 



PCT/FR2004/002104 




RG4 



